Theory of superconductor-ferromagnet point contact spectra: 
the case of strong spin polarization 
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We study the impact of spin-active scattering on Andreev spectra of point contacts between 
superconductors(SC) and strongly spin-polarized ferromagnets(FM) using recently derived boundary 
conditions for the Quasiclassical Theory of Superconductivity. We describe the interface region by 
a microscopic model for the interface scattering matrix. Our model includes both spin-filtering and 
spin-mixing and is non-perturbative in both transmission and spin polarization. We emphasize the 
importance of spin-mixing caused by interface scattering, which has been shown to be crucial for the 
creation of exotic pairing correlations in such structures. We provide estimates for the magnitude 
of this effect in different scenarios and discuss its dependence on various physical parameters. Our 
main finding is that the shape of the interface potential has a tremendous impact on the magnitude 
of the spin-mixing effect. Thus, all previous calculations, being based on delta-function or box- 
shaped interface potentials, underestimate this effect gravely. As a consequence, we find that with 
realistic interface potentials the spin-mixing effect can easily be large enough to cause spin-polarized 
sub-gap Andreev bound states in SC/sFM point contacts. In addition, we show that our theory 
generalizes earlier models based on the Blonder- Tinkham-Klapwijk approach. 

PACS numbers: 72.25.Mk,74.50.+r,73.63.-b,85.25.Cp 



I. INTRODUCTION 

The proximity effect near interfaces between super- 
conductors and ferromagnetic materials has been a field 
of intense research in recent years »i~— This interest is 
mainly triggered by the observation that exotic types of 
pairing symmetries that arc difficult (or impossible) to 
be observed in bulk materials can be created in such 
heterostructures 1 3 i 1 7 i 18 Examples are the recent revival 
of pairing states that exhibit a sign change under the 
exchange of the time coordinates of the particles that 
constitute a Cooper pair ( "odd- frequency pairing"),— or 
mechanisms for the creation of long-range equal-spin 
pairing components in half-metallic ferromagnets.— ~— 
Supercurrents in half-metals have subsequently been 
observed f^ 3 . which ignited a strong activity in further the- 
oretical modeling of this effect &~2il2 r 1 7 1 24 - ^8 gpj n triplet 
pairing has proven to be at the heart of new physical 
phenomena, like 0-7r-transitions in Joscphson junctions 
with FM interlayer o 2 ' 8 ' 12 ' 29 " — or the interplay between 
magnons and triplet pairsi 28 i 32 

So far, transport calculations in SC/FM hybrids have 
mostly been concentrated on either fully polarized FMs, 
so-called half metals (HM), or on the opposite limit of 
weakly polarized systems. However, most FMs have an 
intermediate exchange splitting of the energy bands of 
the order of 0.2-0.8 times the Fermi energy E-p, which 
we here refer to as strongly spin-polarized ferromagnets 
(sFM). As an alternative to solving full Bogoliubov-de 
Gennes equations , 15 ' 26 ^ 3 " — we have recently presented 
a quasiclassical theory appropriate for this intermediate 



range of spin-polarizations, which is of considerable im- 
portance for applications.— 

For such strongly spin-polarized materials, it has been 
argued that Andreev point contact spectra can be used to 
obtain the spin-polarization of the FM,— ~— which is an 
important information for spintronics applications. Ex- 
perimental studies of point contact spectra with strongly 
spin-polarized systems have been performed for a num- 
ber of systems.— ~— However, Xia et a/.— have objected 
rightfully, that without taking into account a realistic de- 
scription of the interface region, the results obtained with 
this method are questionable. 

In the quasiclassical approach to superconducting hy- 
brid structures, interfacial scattering is taken into ac- 
count by the interface scattering matrix S of the struc- 
ture in its normal state. This is ideal for discussing mi- 
croscopic models of interfacial scattering which go well 
beyond the standard Blondcr-Tinkham-Klapwijk (BTK) 
approach.— The latter has been employed to fit experi- 
mental data of SC /FM point contact spectra^ with the 
interface being described by a single parameter Z re- 
lated to its transparency and the ferromagnet by its spin- 
polarization P. The modification of the Andreev point 
contact spectrum compared to a normal metal contact 
is then uniquely related to the spin-dependent density 
of states (DOS) in the FM bulk. This model allows for 
good fits to experimental data. However, comparing dif- 
ferent probes with varying interface transparency, a sys- 
tematic dependence P(Z) was found by Woods et alw^ 
This shows that the extracted spin polarization is not 
a bulk property, as was originally assumed, but at least 
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partially an interface property. This important difference 
has been emphasized also in Ref. l45l 

From the theoretical point of view, it is obvious that 
if scattering is spin-active, i.e. the scattering event is 
sensitive to the spin of the incident electron, this may 
not only imply a spin-dependent transmission probability 
(spin-filtering) (for a review see Ref. |52| ) but also a spin- 
dependent phase shift of the wavefunction^ The latter 
is called the spin-mixing effect and it has been shown to 
be of crucial importance for the creation of exotic pairing 
correlations! 3 ' 19 ' 53 " — So far, estimates of the magnitude 
of this effect and its dependence on physical parameters 
including not only the structure of the interface but also 
the Fermi surface geometry of the adjacent materials and 
the FM exchange splitting are still lacking. Instead, phe- 
nomenological models have been adopted that introduce 
a free parameter to account for it<^ ' 19 ' 54 ' 58 

The main point of this paper is to provide a micro- 
scopic analysis of the characteristic interface parameters. 
In the following we adopt a model of the interface re- 
gion consisting of a spin-dependent scattering potential 
whose quantization axis may be misaligned with that of 
the adjacent FM. We allow for an arbitrary shape of this 
scattering potential and illustrate that this may enhance 
the spin-mixing effect considerably compared to the pre- 
viously used box-shaped or delta-function potentials. We 
also study in detail the relation between spin-mixing an- 
gle and impact angle of the quasiparticlc, showing that 
this relation is non-trivial for transparent interfaces. Fur- 
thermore, we provide a very general mathematical discus- 
sion of suitable parameterizations and representations of 
the scattering matrix in this context. 

Andreev bound states have proven invaluable for 
studying the internal structure of the superconducting 
order parameter i 59 ' 60 Andreev states are also induced at 
spin-polarized interfaces by the spin-mixing effect^ In 
fact, the measurement of such bound states at spin- active 
interfaces would be an elegant method do determine the 
spin-mixing angle of the interface. To date this quantity 
has never been determined in experiment. Our results 
show, that a measurable effect is more likely to appear 
when leakage of spin polarization into the superconduc- 
tor takes place, for example due to diffusion of magnetic 
atoms. Our theory can discriminate between conven- 
tional Andreev reflection processes (AR) and spin-flip 
Andreev reflection (SAR), the latter being responsible 
for the long-range triplet proximity effect. We discuss the 
Andreev bound state associated to the spin-mixing effect 
and show that it may be observable in experiment. Fur- 
thermore we show that for highly polarized FMs, spin-flip 
scattering can influence the spectra considerably, proving 
that such processes must be precluded if one wishes to 
extract the FM spin-polarization from such spectra. 

The paper is organized as follows. In Section [TlJ we 
discuss the quasiclassical theory to describe transport 
through a point contact. In Section Mil we turn to in- 
terface models and discuss the spin-mixing effect and the 
scattering matrix. In Section IIVI we present results for 



Andreev conductance spectra of SC/FM point contacts. 
We dicuss analytical results, focusing on the Andreev 
bound state spectrum, as well as numerical results. In 
Subsection IIV CI we establish the connection to earlier 
transport theories for such systems which are based on 
the BTK approach. We prove analytically that they are 
contained as limiting cases in our formalism. Eventually, 
in Section [V] we conclude on our results. 



II. QUASICLASSICAL THEORY 

We make use of the quasiclassical theory of 
superconductivity 6 -^— to calculate electronic transport 
across the SC/FM interface. This method is based on 
the observation that, in most situations, the supercon- 
ducting state varies on the length scale of the supercon- 
ducting coherence length £o = h\up\/2irk-BT c , with the 
normal state Fermi velocity up. The appropriate many- 
body Green's function for describing the superconducting 
state has been introduced by Gor'kov^ and the Gor'kov 
Green's function can then be decomposed in a fast oscil- 
lating component, varying on the scale of the Fermi wave 
length 2tt /kp , and an envelope function varying on the 
scale of £o- The quasiclassical approximation consists of 
integrating out the fast oscillating component: 



g(p F ,R,e,t) = 



1 



a(p F ) 



&t P hG{p,R,e,t) (1) 



where a(p F ) is the inverse quasiparticlc rcnormaliza- 
tion factor (due to self-energy effects from high-energy 
processes)^ a "check" denotes a matrix in Keldysh- 
Nambu-Gor'kov space,— a "hat" denotes a matrix in 
Nambu-Gor'kov particle-hole space (with 73 the third 
Pauli matrix), p F is the Fermi momentum, R the spatial 
coordinate, e the quasiparticlc energy, t the time, and 
£p = v F ■ {p — Pf)- The quasiclassical Green's function 
obeys the transport equatio n 61 ' 62 



ihvF ■ V^g + [st 3 — A — h, g]t 



0. 



(2) 



Here, A is the superconducting order parameter, h con- 
tains external fields and self-energies due to impurities 
etc, and [•,•]& denotes the commutator with respect to a 
time convolution product (for details see RefHH). Eq. @ 
must be supplemented by a normalization conditio n 61 ' 71 
g ® g = —lir 2 . The current density is related to the 
Keldysh component of the Green's function via: 

j(R,t) = eN F J ^Tr(v F (£ F )t 3 g K (p F ,R,e,t)), (3) 

where N F is the density of states at the Fermi level in 
the normal state, and (•) denotes a Fermi surface average 
which is defined as follows: 



d 2 PF 



N F 



N F J FS (2^) 3 |« F (p F )| 
d 2 PF 

FS (2Trh) 3 \v F (p F )\ ' 



(•), 



(4) 
(5) 
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FIG. 1: (Color online) SC/sFM interface, showing the Fermi surfaces on either side (thick lines). Assuming momentum 
conservation parallel to the interface (k\\), a quasiparticle incident from the SC can either scatter into two (a), or into only one 
(b) spin band of the FM. 



The direct inclusion of an exchange energy Jfm of or- 
der of 0.1 E F or larger in the quasiclassical scheme vi- 
olates the underlying assumptions of quasiclassical the- 
ory. As we aim to describe a strongly spin-polarized FM, 
which means that its exchange field Jfm will be of the or- 
der of the Fermi energy, we cannot include it as a source 
term — ^ Jfm ■ o (with a the vector of Pauli spin matri- 
ces) in the quasiclassical equation of motion. Such an 
approach would neglect terms of order of J FM /E F com- 
pared to A. The resulting condition J F m <^ VE F A, 
assuming e.g. a gap of 1 meV and E-p ~ 1 eV, would 
imply Jfm <C 30 meV. In general, the condition for the 
possibility to include Jfm in the quasiclassical low energy 
scale is violated for most SCs if Jfm > 0.1£f- 

To deal with the strong exchange splitting, we make 
use of the fact that it results in a rapid suppression of su- 
perconducting correlations between quasiparticle states 
with opposite spin, i.e. singlet (|f4-) — |4-t}) ot S x = 
triplet (|t4-) + lit)) correlations. They decay on the short 
length scale Xj = h/(p F2 - Pf3) < hv F2 :i/^ = Co,,- 
Here p F2 , PF3 are the Fermi-momenta of the two spin- 
bands (2 and 3) in the sFM and £or, with r] = 2,3 
the coherence lengths in the respective bands. Conse- 
quently, only equal-spin triplet correlations can penetrate 
the FM-bulk. Hence we pursue the following approach 
to model a strongly polarized FM in the frame of QC 
theory. We define independent QC Green's functions 
(QCGF) for each spin-band which are scalar in spin- 
space, i.e. describe correlations with |yf), respectively 
1 4-4-) spin-wavefunction. The boundary conditions must 
now match three QC propagators at the SC/FM inter- 
face, which we label g v with 77 = 1 = SC, 77 = 2 = f-band 
and 77 = 3 = 4,-band (sec Fig. [TJ. These three QCGFs 
are formally obtained from: 



band separately 



g(p Fv ,R,s,t) 



1 



a(p Fv ) 



d£ pv T 3 G(p,R,e,t) (6) 



with £ pn = vp v ■ (p — p Fv ), PFr] and vf v being the re- 
spective Fermi-momenta/velocities of the bands. Con- 
sequently, the current must then be evaluated for each 



de 



j v (R,t) = eN Fn J — Tr^VFniPF^hgr, (PF V ,R,£,t) 

(7) 

Here, Nf v is the partial density of states at the Fermi 
level in band 77, and (•) r; denotes the corresponding Fermi 
surface average 



(•>, 



1 



d 2 pFr, 



Nf v Jfs v (2Trh) 3 \vFn(PFr,)\ 

d 2 pF V 

FS V (2TTh) 3 \v Fll (p Fr ,)\ ' 



(•) 



(8) 
(9) 



In addition, the system's properties vary on the atomic 
length scale in the interface region between the two ma- 
terials. Thus the QC theory is also not applicable in 
the immediate proximity to the interface (on the scale 
of the Fermi wavelength). This is a general problem in 
the quasiclassical description of heterostructures, which 
can be circumvented by deriving appropriate boundary 
conditions for matching the QC propagators on both 
sides of the interfaced The full boundary conditions for 
the present problem have been developed only recently^ 
Earlier works on Andreev spectra using QC theory were 
restraint to either SC/normal metal contacts with spin- 
active interfaces , 45 i 54 ' 73 i 74 or contacts with weak ferro- 
magncts. We refer to Refd and references therein for a 
detailed discussion of this problem. In the following sub- 
section we discuss a parameterization of the QC propa- 
gator, and return to the problem of boundary conditions 
at the interface in Subsection IHBI 



A. Riccatti parameterization 

For our calculations we choose a representation of the 
quasiclassical Green's function (QCGF) that has proven 
very useful in the past and is standard by now. In this 
representation, the Kcldysh QCGF is determined by six 

of 



parameters in particle-hole space, 7 



R,A ~R,A „K ;K 

5 / 
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which 7 R,a (pf, R, £, t) and 7 R,a (pf, R, £, t) arc the re- 
tarded (R) and advanced (A) coherence functions, de- 
scribing the coherence between particle-like and hole-like 
states, whereas x K (pF, R, e, t) and x K (pp, R, e, t) are dis- 
tribution functions, describing the occupation of quasi- 
particle states.—^ The coherence functions are a gener- 
alization of the so-called Riccatti amplitude s) 77 ' 78 to non- 
equilibrium situations. All six parameters are 2x2 spin- 
matrix functions of Fermi momentum, position, energy, 
and time. The parameterization is simplified by the fact 
that, due to symmetry relations, only two functions of 
the six are independent. The particle-hole symmetry is 
expressed by the operation (•), which is defined for any 
function of the phase space variables by 

Q(p F ,R, z, t) = Q(-0r, R, -z*,t)*, (10) 

where z — e is real for the Keldysh components and 
z is situated in the upper (lower) complex energy half 
plane for retarded (advanced) quantities. Furthermore, 
the symmetry relations 

7 A = (7 R ) f , 7 A = (7 R ) f , * K = (* K )t (11) 

hold. As a consequence, it suffices to determine fully the 
parameters 7 R and x K . 

The QCGF is related to these amplitudes in the fol- 
lowing way [here the upper (lower) sign corresponds to 
retarded (advanced)] J£ 

g R ' A = T2m ( | ±MTf 3) (12) 

with the abbreviations Q = (1 — 77) _1 and T = Qj, and 

^-K-UT(o K J0(-U) A - 

(13) 

Note that all multiplication and inversion operations in- 
clude 2x2 matrix algebra (and, more general, for time- 
dependent cases also a time convolution). 

From the transport equation for the quasiclassical 
Green's functions one obtains a set of 2x2 matrix equa- 
tions of motion for the six parameters above i 75 ' 79 For 
the coherence amplitudes this leads to Riccatti differen- 
tial equations* 7 - 8 - hence the name Riccatti parameteriza- 
tion. As we are interested in this paper only in the inter- 
face problem in relation to a point contact, the transport 
equations are not relevant for the problem at hand. For a 
point contact, the superconductivity is modified only in a 
very small spatial region, and this modification can be ne- 
glected consistent with quasiclassical approximation. We 
assume that the half-space problem is solved and calcu- 
late the conductance across the point contact. For this, 
we turn now to the problem of solving the boundary con- 
ditions for the point contact. 



B. Boundary conditions 

1. General case 

The QCGF mixes particlelike and holelike amplitudes, 
and as a result the transport equations are numerically 
stiff, with exponentially growing solutions in both posi- 
tive and negative directions along each trajectory, which 
must be projected out. A particular advantage of the co- 
herence and distribution functions is that, in contrast to 
the QCGF, they have a stable integration direction for 
each trajectory. This direction coincides with their prop- 
agation direction, and is opposite for holelike and parti- 
clelike amplitudes as well as advanced and retarded ones. 
This allows to distinguish between incoming and outgo- 
ing amplitudes at the interface. We adopt the notation 7 -^ 
that incoming amplitudes are denoted by small case let- 
ters and outgoing ones by capital case letters. Boundary 
conditions express outgoing amplitudes as functions of 
incoming ones and as functions of the parameters of the 
normal-state scattering matrix. They are formulated in 
terms of the solution of the equation^ 

[iw] R = b' kk . + £ n^ kl %^' klk ,] R (14) 

for [rfc.j-fc/] 11 , where the trajectory indices k,k',ki run 
over outgoing trajectories involved in the interface scat- 
tering process, and the scattering matrix parameters en- 
ter only via the "elementary scattering event" 

p 

(the trajectory index p runs over all incoming trajec- 
tories). It is useful to split the quantity [Ffc^fc/] R into 
its forward scattering contribution, which determines the 
quasiclassical coherence amplitude, 

r R = r^ fc , (16) 

and the remaining part 

[IV fc ,] R = [iw -r fe 4fc<] R , (17) 

which is relevant only for the Keldysh components. Anal- 
ogous equations^ hold for the advanced and particle- 
hole conjugated components, [r p ^ p /] R , [r p '_> p ] A , and 
[r\.'^fc] A . The boundary conditions for the distribution 
functions reacU^ 

xf = £ [6 kkl +r^ fel7fcl ] R [4 ifc2 ] K [ ( 5 fc2fc + 7fc2 r w ] A 

fci ,k 2 

-^[r^ fel ] R i R [r fcl ^] A , (is) 
fci 

which depend on the scattering matrix parameters only 
via the elementary scattering event 

[ x fcfc'] K = £ -Sfcp^p Sp k '- (19) 
p 

Analogous relations hold for A R . 
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2. Special case: SC/FM point contact 

In the case under consideration the trajectory labels k 
and p run from 1 to 3, with 1 denoting (spin-dcgcncrate) 
trajectories on the superconducting side, and 2 and 3 tra- 
jectories for the two spin directions on the ferromagnetic 
side. We use the following notation for the (unitary) 
scattering matrix: 





(Ri 


Tl2 


T 13 \ 


*- 




T 2 






\T 31 


r-ii 


r 3 J 



The current across the interface is conserved (this is en- 
sured by our boundary conditions), so that it suffices 
to calculate the current density at the FM side of the 
interface. We proceed with expressing the outgoing am- 
plitudes for bands 2 and 3 in terms of the incoming am- 
plitudes and the scattering matrix. 

For a point contact with semi-infinite SC and FM 
regions (assuming that the Thouless energy related to 
the geometry of the system is negligibly small), there 
are no incoming correlation functions from the FM side, 
7^3 A = 723 = 0, whereas on the SC side we can use the 
bulk solutions. For a singlet order parameter the bulk 
solutions of the coherence functions read 

R,A A s Zff 2 ~R,A _ A*Z<7 2 

71 e±VIA s | 2 -e 2 ' 71 " e±VI A sl 2 -e 2 ' 

(21) 

with the singlet superconducting order parameter A s . 
Taking into account these facts, we obtain from Eq. (fT4)) 

T^i = [72i] R + r 2 ti7i R [7ii] R , (22) 
= [72 2 ] R + r 2 ti7i R [7i 2 ] R ) (23) 

r^a = [7 23 ] R + r 2 Vi7 R [7i 3 ] R , (24) 

with [ 7 y R = Sajf-SZj for i,j = 1,2,3. The first equa- 
tion, Eq. (|22]l. can be solved, 

r£_i = T 2ll ?R{ (1 - jf R^fRiy 1 . (25) 

It appears useful to introduce the notation 

A = Tf^f. (26) 

From Eqs. (|23 ]) -(|23 jl we obtain 

Tf = (T 21 + AR 1 )^T* 2 , (27) 
r£- 3 = (T 21 + AR.hfT*,. (28) 

Note that the identity T 2 i+AR 1 = T 21 (l--ffR^R 1 )- 1 
holds. The corresponding solutions for band 3 are sim- 
ply obtained by replacing 2 O 3. Amplitudes and 
f£_ 3 are obtained using Eq. ([TO]), with = S*j. The 
required advanced amplitudes can be obtained from the 
fundamental symmetry relations of this formalism, which 
imply f£ = (r*)t and f^ a = (T^ . 



For the distribution functions, we use a gauge in terms 
of anomalous components.— Taking the electrochemical 
potential equal to zero in the SC, and equal to — cV in 
the ferromagnet, these are x\ = x\ =0 and 



^2,3 = 
X2,3 = - 




Note that in our notation e = — |e|. From Eq. (fTS)) we ar- 
rive at the following expressions for the outgoing Kcldysh 
amplitudes for band 2: 

x 2 = [4 2 ] K + r 2 t 1 7i R K 2 ] K + [4i] K 7^ 2 
+r 2 Vi7iVn] K 7i i f t* - r^ 3 ^f ^2 (30) 

with [x' tj ] K = S i2 x 2 Sj 2 + S i3 x 3 Sj 3 for i,j = 1,2,3. In- 
troducing what has been obtained before, we arrive at 

X 2 = {r 2 +AT 12 )x 2 {r 2 +AT 12 y 
+ (^23 + AT 13 )x 3 {r 23 + AT 13 ^ 
-(T ai + AR 1 )( 7 fT^ 3 )x 3 ( 7 fT^(T 21 + ARtf (31) 

Again, the corresponding solution for band 3 is obtained 
by replacing 2«3. 

III. INTERFACE MODEL 

We consider a point contact with a diameter much 
smaller than the superconducting coherence length but 
still larger than the Fermi- wavelength, as shown in Fig. [2] 
a. A larger contact would result in a perturbation of 
the SC state, a smaller one would invoke conductance 



(a) (b) 




FIG. 2: (Color online) (a) The Andreev point contact with 
spin-active interface; (b) The spherical angles a and ip char- 
acterize the orientation of the interface exchange field Ji with 
respect to the FM exchange field Jfm- The dashed arrow 
indicates the area where the misaligned interface magnetic 
moment resides. 
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quantization^ -. This also allows for the decisive assump- 
tion of translational invariance on the scale of the Fermi 
wavelength Ap. The region in the immediate vicinity of 
the interface (I) cannot be described within QC theory. 
Instead, the normal state scattering matrix of the inter- 
face must be obtained from microscopic calculations and 
then enters the QC theory through boundary conditions 
as outlined above. 

The mechanism giving rise to spin-active scattering at 
the interface is the ferromagnetic exchange held in both 
the adjacent ferromagnetic material and in the interface 
itself. The interface will in general carry a magnetic mo- 
ment, that in the simplest case is induced by the mag- 
netization of the bulk ferromagnetic material; however, 
there might be cases where an extra interface magnetic 
moment develops, either manufactured by using a thin 
magnetic layer, or due to spin-orbit coupling, and related 
to that, magnetic anisotropy. The interface magnetic mo- 
ment can be misaligned with the one of the bulk sFM. We 
characterize this misalignment by two spherical angles a 
and ip, as indicated in Fig. [5] b. While the spin-activity 
of interfaces has been discussed extensively in the theory 
of superconducting heterostructures, most of the work 
so far considered a set of phenomenological parameters 
for characterizing the intcrfacial scattering. Notably, one 
of these parameters, the so-called spin-mixing angle, or 
spin-dependent phase shift, turned out to be of decisive 
importance for the creation of unconventional supercon- 
ducting correlations in proximity to the interface. The 
spin-mixing angle is essentially a relative phase differ- 
ence between f and I electrons acquired upon scattering. 
Obviously, an exchange field in the interface region will 
provide such an effect, but other mechanisms, like for 
instance spin-orbit coupling arc also candidates. 

So far, estimates of the possible magnitude of this ef- 
fect based on a physical model of the interface region 
are still lacking. Here, we will provide such an analysis 
based on wavefunction matching techniques. In partic- 
ular, we will discuss the dependence of the spin-mixing 
effect on the shape of the barrier. To this end, we con- 
sider a spin-split potential barrier which is assumed to 
conserve the momentum component parallel to the inter- 
face upon scattering. For the system we deal with, this 
gives rise to two types of transmission events (see Fig. [TJ . 
Depending on the impact angle the parallel momentum 
conservation constraint either allows for or prohibits scat- 
tering that involves the minority spin-band of the sFM. 
For a half metal, where the 4,-band is completely insulat- 
ing, only the latter case occurs. 



A. Interface scattering matrix 



all of these parameters will be relevant for the physical 
problem at hand. For instance, spin-scalar phase factors 
do only matter for two or more interfaces. Furthermore, 
since a singlet SC is spin-isotropic, one is free to choose 
the spin-quantization axis in the SC conveniently. To 
clearly identify these irrelevant parameters we use a spe- 
cial parameterization of a general unitary matrix with 
the aforementioned dimensions, as discussed in the Ap- 
pendix. The most important result of these considera- 
tions is that the spin-mixing effect can be fully described 
by only one parameter in the HM case, but three are 
required in the FM case. 

Neglecting irrelevant spin-scalar phases and using the 
gauge freedom in the SC the scattering matrix reads for 
the first type of scattering 



S 



FM 



ii5/2 



t 2 e 



t' e 



t 2 e li} ^ 2 
t' 2 e-^/ 2 



'2 



r 2 3 



r23 r 3 j 

(32) 

The scattering matrix for the second, HM type, scatter- 
ing is 



(33) 









t 2 e^/ A 


Shm — 




r n e- M ' 2 


t' 2 e~ M / 4 




^ t 2 e^ 




r 2 



There is also the possibility of total reflection with no 
transmission on either side, in which case the scattering 
matrix consists of the reflection parts only. In writing 
the scattering matrices (|3"2"|) and (|3"5")) we have put the 
(p-phasc that appears in Fig.[2][b) to zero, since the prob- 
lem we consider is invariant with respect to rotation of 
the interface magnetic moment around the bulk magne- 
tization; the scattering matrix is symmetric in this case, 
S = S T . We also omitted possible complex phases in 
the reflection part on the FM-side, i.e. r 2 , and r 2 ^, as 
they are irrelevant to the problem at hand. The require- 
ment of unitarity leads to additional relations between 
the reflection and transmission parameters. The phases 
that we wrote explicitly in Eqs. (|3"2")) and ([33)) are cru- 
cial, since they account for the spin-mixing effect. In the 
following section, we will discuss their magnitude as a 
function of various interface parameters. 

Using the set of independent parameters described in 
the appendix we have: 



ri t = r t cos(a y /2) 2 + r x sin(a y /2) 2 
r u = r t sin(ay/2) 2 + 7-4 cos(a y /2) 2 
sin(ay) 

r m = -( r t-n) — ~ — • 



(34) 



At this point we present some general considerations 
concerning the scattering matrix of a spin-active inter- 
face. Such a matrix is unitary and of dimensions 4 x 4 in 
the FM and 3 x 3 in the HM case. The maximum number 
of free parameters is 16 or 9 respectively. However, not 



The angle ay defines a rotation in spin-space to the inter- 
face eigenstates, characterized by transmission and reflec- 
tion eigenvalues. Its precise definition is given in the ap- 
pendix. Most importantly, it is in general not identical to 
the interface misalignment angle a, however approaches 
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it in the limit of thick interfaces. For thin interfaces it 
is renormalized by the influence of the exchange field of 
the adjacent FM. r-j- and are the singular values of the 
reflection block R\. In the tunneling limit, r-j-, rj. w 1, 
and the off-diagonal elements vanish even for ay 7^ 0. 
This is easily understood from a physical point of view, 
since spin-flip reflections on the SC side require that the 
reflected quasiparticles "feel" both misaligned exchange 
fields, J\ and Jfm, and not just that of the interface. It 
is possible to provide analogous expressions for the re- 
maining parameters of the scattering matrix, however in 
the sFM case they are rather cumbersome and also not 
needed for the following analytical discussion. For the 
half-metallic case, the only relevant phase parameter is 
the spin-mixing-anglc #, and for the remaining parame- 
ters we have r\. = 1 and 



t-f cos 



(t)-« 



tf sin 



(^) , r 2 = -r t . (35) 



In the following we will discuss the influence of the 
shape of the scattering potential, and will show that the 
widely used box shaped or delta-function shaped poten- 
tials gravely underestimate the magnitude of the spin- 
mixing effect. 



B. Box-shaped scattering potential 

In this section we consider spin-dependent box poten- 
tials, for which analytical solutions can be obtained. In 
particular, we discuss here the dependence of the spin- 
mixing angle, on the impact angle of the incoming 
quasiparticle which is parameterized by the momentum 
component parallel to the interface, kn. The model pa- 
rameters are the misalignment angle a (sec Fig.[5]b), the 
energies of the band minima in the FM with respect to 
that in the SC (£2, £3), the spin-dependent height of 
the potential (U+, U-), and the width of the potential d 
(see Fig. [3]). All energies are given in units of £f and d 
in units of Af/27t with Ap = 27r/fcp the Fermi wavelength 
in the superconductor. 

The scattering matrix is defined with respect to the 
chosen spin-quantization axes on both sides of the inter- 
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FIG. 3: (Color online) Sketch of the box-potential model that 
we consider in this section (right) and of the corresponding 
Fermi-surface geometry (left). The model parameters are in- 
dicated. 



face. Naturally, on the FM side we use the bulk sFM 
magnetization axis. On the SC side we use that of the 
interface magnetic moment. To obtain an S-matrix with 
the structure defined above, one must subsequently cal- 
culate and apply a rotation of the quantization axis in 
the SC: 







s 



Q 
1 



(36) 



where Q is a spin rotation matrix acting on spins in the 
superconductor. We describe this procedure in App. [A~T1 
All the quantities plotted are calculated in this rotated 
frame, the point being that otherwise one does not have 
an unambiguous definition of the mixing-phases. Nat- 
urally, the Andreev spectra are invariant under these 
transformations. We obtain the scattering matrix by 
matching wave functions as described in App. IA 21 

In Fig.|3]a,b we show the spin- mixing angle for different 
values of the interface potential width d. The band min- 
ima in the FM are E 2 = 0.1 E F and E 3 = 0.9 E F , which 
implies that at k\\ > 0.31 /cfi the minority band becomes 
insulating and the scattering matrix reduces to a 3 x 3 
matrix. In the tunneling limit (d 3> Af/27t) the spin- 
mixing angle behaves as expected: it is approximately 
given by the value (see App. IA 2|) 



d = 2 



which appoaches zero for grazing impact (fci 



(37) 
0), and 



2[arctan v /£>/([/+ - £' F )-arctan A /£ , F /(t/_ - E F )] (« 
0.297T for Fig. [4} for normal impact. Here, fei is the com- 
ponent of the wavevector perpendicular to the interface 
in the superconductor, and k± are the exponential decay 
factors for the spin- up/down wave functions in the bar- 
rier. For thin (highly transparent) interfaces the mixing- 
angle $ is a more complicated function of the quasipar- 
ticle impact angle. In this regime, 1? is predominantly 
controlled by the Fermi-surface geometry indicated in 
Fig. [3J There is a local minimum at k\\ > kp 3 , and for 
very thin interfaces $ is largely enhanced for grazing im- 
pact (d = 0.1 Ap/27r in Fig. 0]). This enhancement can 




FIG. 4: (Color online) The spin-mixing angle 1? as function 
of the momentum component parallel to the interface, shown 
for various barrier thicknesses, (a) d — 0.1, 0.5, 1.0 Af/27t, 
(b) d — 2.0, 3.0, 5.0 Xf/2tt. The remaining parameters are 
E 2 = 0.1.Ep, E z = 0.9 E F , U+ = 1.1 E F , U- = 1.9 E F , 
a = 0.5 7r (see text and Fig. O- 
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be understood from the d — limit, i.e. the case where 
the interface barrier is absent. In this case, 



d 



2 arctan 



(38) 



where K3 corresponds to the imaginary wave vector in the 
insulating band 3, which controls the exponential decay 
of the spin-down wave function into the fcrromagnet. In 
the particular case we show here, see Fig. [3J ki takes 
a finite value for all trajectories that contribute to the 
current, while K3 increases monotonously from at fen = 
fcp3 ~ 0.31 fcpi to some finite value at k\\ = fcp2- This is 
because the effective height of the potential for tunneling 
into the insulating band increases with k\ \ . For Fcrmi- 
surface geometries with fcpi < fcp2 (not shown here) the 
wave vector k\ drops to zero for grazing impact, and so 
does the spin-mixing angle. 

In the present case, the situation is complicated by 
the fact that we consider both a finite interlayer and a 
broken spin-rotation symmetry. This leads to a finite 
spin- mixing angle even for k\\ = &F3 and below, which 
leads to the non-trivial behavior with a minimum for in- 
termediate impact angles. This illustrates that not only 
the scattering potential itself but also the Fermi-surfacc 
geometry is highly important for spin-active scattering 
beyond the tunneling limit. 

As for the magnitude of the mixing effect, we stress 
that for a realistic choice of parameters, it is hardly pos- 
sible to achieve mixing-phases above 0.57T in this model. 
In Fig.@]we use an exchange field of J = 0.8 Ep, which is 
close to the half-metallic limit. Using smaller exchange 
energies naturally leads to a smaller effect, as can be seen 
in Fig[5] a,b, where we plot "& for different values of the 
exchange field J = E3 — E2. 

In Fig. [5] we show the spin- mixing phases associated to 
transmission, $2 and $3. One can see that $2 = #/2 for 
fen > 0.31 fcFi- This relation one would expect for a SC 
contacted with a half-metallic ferromagnet; the finding in 
Fig. [5] is consistent with this and the discussion presented 
above, since the trajectories under consideration effec- 
tively correspond to the HM case. For ku < 0.31 fcpi, the 
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FIG. 6: (Color online) The spin-mixing angles #2(fc||) [(a) 
and (b)], and #3 (fen) [(c) and (d)], for thin (left column: d = 
0.1 (solid), 0.5 (dashed - dotted), 1.0 (dashed) A f /2tt) and 
thick (right column: d — 2 (dashed), 3 (dashed — dotted), 
5 (solid) Af/2-7t) interfaces. The remaining parameters are 
the same as in Fig. [3] 



mixing phase is considerably enhanced above the value 
oi'd/2. The plots also illustrate that i?2 and #3 are differ- 
ent in magnitude and also vary differently with ku . As we 
show in the appendix, the mixing-phases $2 and $3 are 
correlated with $ but in general also depend on a number 
of other free parameters. Their magnitude is decisive for 
the creation of triplet correlations in the corresponding 
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FIG. 5: (Color online) (a) Spin-mixing angle ■& as a function 
of impact angle for (a) d = 0.5 Xf/tt, and (b) d — 5.0 Af/27t. 
In both plots, the curves are for U- = 1.2. ..2.0 Ef, Es — 
U- — 1.0. The corresponding value of the exchange field J is 
indicated. The remaining parameters are E% = O.ISf, U+ — 
1.1 Ef, and a = 0.5 7T. 



FIG. 7: (Color online) The transmission parameters |t2£2|(k||) 
[(a) and (b)], and |i3t 3 |(fc||) [(c) and (d)], for thin (left column: 
d = 0.1 (solid), 0.5 (dashed - dotted), 1.0 (dashed) Af/2tt) 
and thick (right column: d = 2 (dashed), 3 (dashed — dotted), 
5 (solid) A F /27r) interfaces. The remaining parameters are the 
same as in Fig. f3] 
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band, as we will show below. 

In Fig. [7] we present the product \t v t'A (which con- 
trols the magnitude of long-range SAR). We plot this 
quantity for both the majority (upper row) and minority 
(lower row) band of the FM. Apparently there is a non- 
monotonous dependence on the interface width d, which 
is related to the fact that spin-flip scattering becomes 
more effective as the interface region becomes larger. For 
even larger d the global suppression of transmission in- 
tervenes and we approach the tunneling limit. Again, we 
note that for thin interfaces the dependence on trajec- 
tory impact angle is non-monotonous, showing maxima 
for non-perpendicular impact. These maxima coincide 
exactly with the minima of the spin-mixing angle. Note 
that a nonzero tj ; requires a non- vanishing misalignment 
angle a. 

To conclude on this section, we have shown that the 
magnitude of the spin-mixing effect is limited to rather 
small values in the box potential case if one assumes 
J < E-p and d s Ap. Moreover, both spin- mixing ef- 
fect and spin- flip scattering are very sensitive to trajec- 
tory impact, interface thickness, exchange field of the 
interface and the Fermi surface geometry of the adjacent 
materials. 



C. Delta- function scattering potential 

A special case of the box-shaped potential is that of the 
delta-function potential, that is widely used in describing 
interfaces within the BTK paradigm. Here, we show that 
the situation is in this case comparable to that of the 
box potential. Delta-function models introduce a weight 
factor Vq of the Delta- function which enters the matching 
condition for wavefunction derivatives: 

±V 1 {z = 0)-±* 2 {z = 0) = ?^9 a (z = 0). (39) 

A spin-dependent potential can simply be modeled by 
choosing a spin-dependent weight factor V± . This weight 
factor effectively corresponds to the area under the scat- 
tering potential, i.e. we have V± = (U± — Ep) ■ d, to 
connect with the notation above. In Fig. [5] we plot # as a 




0.00 0.25 0.50 0.75 1.00 0.00 0.25 0.50 0.75 1.00 
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FIG. 8: (Color online) (a) Spin-mixing angle $ as a function 
of V+/V- for the Delta- function potential. E2 = 0.1 Ef, 
E z = 0.7 E F . (b) The same as (a) for E 2 = -0.7 E F , E 3 = 
-0.1 E F . 



function of V+ /V— for perpendicular impact and two dif- 
ferent choices of the Fermi-surface geometry. Since we do 
not calculate any spectra for this model, we choose a = 
for simplicity. Generically, spin- mixing angles d > 0.57T 
can only be reached for V+/V- < 0.1, which requires ei- 
ther V+ to be very small or an interface exchange field 
exceeding the Fermi-encrgy. 



D. Scattering potentials with arbitrary shape 

The box potential actually constitutes a high degree of 
idealization. The most obvious generalization is to con- 
sider a potential that varies smoothly on the scale of a 
few interatomic distances, or on the scale of the Fermi 
wavelength in metals^ This is quite realistic taking into 
account that metallic screening of charges takes place 
only on the Thomas-Fermi wavelength scale. In addition, 
some magnetic ions might penetrate the superconductor 
from the ferromagnet, leading to a spin-dependent po- 
tential that decays in the bulk of the superconductor. 
In the latter case a certain degree of disorder is intro- 
duced. However, we will assume that any such disorder 
is weak, so that the momentum component parallel to 
the interface is still a good quantum number. A truly 
realistic description would have to drop the assumption 
of translational invariance and consider disorder on a mi- 
croscopic level. In principle our theory can be extended 
to this regime, but this is beyond the scope of this paper. 
If one is only interested in transmission and reflection 
amplitudes, the difference between the box-shape and a 
smoothened potential is negligible. But when scattering 
phases are important, as in our case, this is not true, as 
we will show in the following. 

For definitcness, we consider a potential shape as 
shown in Fig. [5J with Gaussian "slopes". The "smooth- 
ness" of the interface barrier is then controlled by the 
standard deviation a of the Gaussian. Hence, we have 
the spin-dependent potential: 

f (Vi ± J/2) • e -(2+0> 2 z <-d 
U± = < Vi ± J/2 -d < z < . 

[ E n + (Vi ± J/2 - E ri )e~ z2 ^ 2 z > 

(40) 

In the limit of a very smooth potential, one may resort to 
the Wentzel-Kramers-Brillouin (WKB) approximation^ 
to calculate the scattering problem. An interface that 
complies to the requirements of WKB would have to be 
much larger than the Fermi-wavelength however, which 
is unrealistic. For this reason we resort to a numerical 
method for calculating the scattering problem. We use 
a recursive Green's function technique^ to calculate the 
single particle Green's function of the interface Hamilto- 
nian and obtain the scattering matrix from it using the 
Fisher-Lee relations.— To study the effect of the poten- 
tial shape on the spin-mixing angle, we plot the angle i9 
in Fig. [TO] b for different values of a. To avoid a large 
variation of the interface transmission when varying er, 
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FIG. 9: (Color online) Sketch of the scattering potential for 
the smooth potential model (right) and the here considered 
Fermi-surface geometry (left). The parameters introduced in 
Eq. 2n]are indicated. 



return points for the respective spin bands (see Fig.[3]for 
the notation). In the intermediate case that we consider 
here both the different wavevector mismatches and the 
dephasing of propagating modes will add to the mixing 
effect. Note that only the shape on the SC side of the 
interface matters for boosting the spin-mixing angle 1?. 

The discussion in terms of scattering matrix parame- 
ters presented here is flexible enough to be extended, e.g. 
to other Fermi surface geometries, or adiabatic variation 
of the interface magnetization. Furthermore, instead of 
insulating interfaces one could consider interfaces where 
one or even both channels are conducting. The latter 
case has been considered by Beri et alJ^ 



we keep d + a = 0.7 Af (see Fig. HUla). 

Furthermore, we use E 2 , E 3 < here, i.e. both the 
FM-bands have a larger Fermi-surface than the SC. As 
we will see later on, this Fermi surface geometry and the 
scattering constraints it implies can have an important 
effect on the shape of the spectra, and in particular on 
features which arc related to the spin- mixing effect. 

The main result of considering a variation of the po- 
tential shape is however, that it has a tremendous effect 
on the spin- mixing angle, as clearly seen in Fig.[TU]b. Its 
magnitude can exceed for a smooth potential that for a 
box potential of similar transmission easily by a factor 
of 3-4 or more. This is sufficient to observe some ex- 
otic features related to this effect in the Andreev spectra 
of point contacts, as discussed in the next section. The 
physical reason for this is that, unlike in the box potential 
case, electrons with opposite spins acquire a phase differ- 
ence while they are still propagating, which implies that 
a larger mixing phase is not inevitably tied to a strongly 
reduced transmission. This can be best seen in the WKB 
limit, where the mixing angle is exclusively given by this 
dephasing: 



= 2 



dzp^{z) 



(41) 



Here p-^j, = ^/2m(Ep — U±) and Zf t ± are the classical 




FIG. 10: (Color online) (a) Shape function of the scat- 
tering potential (average between both spin directions) for 
a = ... 0.7 A F and a + d = 0.7 Af. E 2 = -0.1 E F and 
E3 = —0.8 Ef. (b) The spin-mixing angle $ as a function 
of impact angle for the different potentials plotted in (a), a 
increases from bottom to top. 



IV. ANDREEV CONDUCTANCE SPECTRA OF 
SC/FM POINT CONTACTS 

In the remaining part of the paper we discuss Andreev 
spectra that result from our model. We use a definition 
for the FM's spin-polarization given by 



P 



N F2 - N F3 



F2 



1-':! 



(42) 



For parabolic bands, the density of states is proportional 
to the Fermi-momentum, N Fri oc p Fv cx y/Ep — E v . We 
assume equal effective masses. 

The current density in terms of the distribution func- 
tions and coherence functions is given by 



Jn = 



Frj 



de (v v ■ jsj-i) v +, 



(43) 
(44) 



where the expression for j e>ri is given by 



Je:2 = x 
—x 



{|r 2 + AT 12 \ 2 + |r 23 + AT 13 \ 2 - l}- (45) 
{|(T ai + ARrX^TTg)! 2 + |(T 21 +^ J R 1 )( 7l T* 2 )| 2 }, 



and an analogous expression is obtained for j e 3 by in- 
terchanging 2 with 3. Here, (»), ; + means a Fermi-surface 
average over one half of the Fermi surface (positive mo- 
mentum directions, pointing into the FM, for the first and 
third term of Eq. (|4"4")) . negative directions for the second 
term). To derive this expression, we used the universal 
symmetry relation (|10[) . Furthermore x = x 2 = x 3 as 
defined in (|2U)) . A is defined in (f2l))) and the scattering 
matrix parameters in (|2"0")) . Equations (|4"3"|) -([4"5 j) are the 
main result of this paper. 

The interpretation of Eqs. (|3T|) and f|44[) allows for 
identifying two types of Andreev reflection, shown in 
Fig. II H one of them giving rise to a long-range prox- 
imity effect in the FM. The terms T^x 2 f £ in Eq. (gj) 
and r 2< _ 3 x 3 T 3 ^ 2 entering X 2 in Eq. ([30]) both describe 
current contributions from Andreev reflected holes to the 
current in band 2. The first term is proportional to the 
incoming distribution function in the same band. Thus 
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(b) 



(d) 





FIG. 11: (Color online) Transport processes contributing to 
the current through the point-contact (a) Normal Andreev 
reflection (AR) (b) Normal transmission, requires e > A (c) 
Spin-flip Andreev reflection (SAR) (d) Spin-flip transmission, 
requires e > A 



we refer to it as spin- flip Andreev reflection (SAR), as it 
requires a spin-flip to transmit a singlet pair into the SC. 
The second term corresponds to normal Andreev reflec- 
tion, since it reflects a hole in the opposite band. While 
SAR is related to the outgoing (equal-spin) triplet corre- 
lation function in the respective band, AR is described as 
a term renormalizing the outgoing distribution function. 
Unlike SAR, AR does not contribute to the coherence 
functions in the FM spin-bands. 

Using the scattering matrix parameterization intro- 
duced in Sec. IIII1 we can obtain explicit analytical so- 
lutions for the coherence functions: 



r R = 

V 



T vl Nia y T? v 
D 



T 3 iNia y T 12 



D 



(46) 



= 7 (1 + 77 [e-***' (p + Ps) - Ps + 2 P a y sin(0/2)] ) 



D =1 + 2jj[pcos$ — p s (l — cos?9)] + (77) 2 p 2 



P = r t r l Ps = ( r T - r l) 



2 sin 2 (ay) 



(47) 
(48) 

(49) 



p = v / pl • diag[r-f cos 2 (ay/2) + sin 2 (ay/2), 
r± cos 2 (ay /2) + sin 2 (ay /2)] 

Here r\ £ {2, 3}, and we omitted the index 1 for the in- 
coming coherence functions. -f R is related to 7^ in (f2~Tj) 
by j R i<7 y = J R . The advanced component is obtained 
via f r f = (!*)•. 

Note that the T-functions differ only by the transmis- 
sion vectors T, but since the numerator is a matrix prod- 
uct, this still gives expressions that differ markedly. In 
any case, we have = if ay = or if •& and i? I( = 0. 
We focus on the denominator, Eq. (|48|). which arises from 
the matrix inversion in Eq. (|25[) and is the same for all 



coherence functions. It is of particular interest since it 
leads to the emergence of conductance peaks in the An- 
dreev spectrum. 



A. Andreev bound state spectrum 

The appearance of the Andreev conductance peaks can 
be seen most clearly in the tunneling limit. Here p s = 
and p = 1 which simplifies the expressions above consid- 
erably. The full solutions read 



" " 7 l + (77) 2 + 2cos(tf)77 



2<-3 



t,j77 



1 + (77) 2 + 2 cos(i?)77 



(50) 
(51) 



with 



r R =t 2 t 3 e^ 



'•23 t 32t 



I23t32e 3 



^23 = (l?2 + 1? 3 )/2. 



For e < A we have 7 



"7 



and 1 7 fl | = 1 and we 
can easily show that (|50|) and ([5"Tj) both have a pole at^ 



(52) 



e po ic = ±Acos(tf/2). 



This pole corresponds to an Andreev bound state induced 
by the spin-mixing effect at the superconducting side of 
the sample. Following Fogclstrom,— one can show that 
these bound states appear in the DOS of the supercon- 
ductor close to the interface and are actually associated 
to different spins. For tt > 1? > 0, the bound state for 
e > appears in the DOS of ^-quasiparticles and that 
for e < in that of J,-quasiparticles. The appearance 
of the sub-gap peak is only tied to the spin-mixing an- 
gle It does not depend on spin- flip scattering or the 
mixing phases associated to transmission. However, we 
shall see that a high mixing angle of > 0.5 n is required 
to make this bound state appear in a finite temperature 
spectrum. If we consider the full expression of the de- 
nominator, Eq. (|48[) . we find that the pole is lifted, yet 
2 local maxima remain (see Fig. I12[) which decrease in 
magnitude with increasing transparency of the interface, 
as the bound state acquires a finite lifetime. Obviously, a 
tunneling barrier in addition to an appreciable mixing ef- 
fect is required to observe a sub-gap-peak in the Andreev 
spectrum. 

In the half-metallic case, the full solution for the out- 
going coherence function reads: 



-?Tat 2 (l + r)7 R (l~7 R 7 R ) 
1 + 7 R 7 R[2r cost? - T^t 4 } + (jRjR^r 2 



(53) 



with Ta = sin($/2) sinfay )/ (1 + r). This solution is al- 
ready discussed in Ref. [la, but we state it here again 
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FIG. 12: The denominator of equation (|46|l as a function of 
the quasiparticle energy e. r-f = 0.9, rj_ = 0.95, ay = 0.5 tt. 
Plots for i? = 0.1 7T...1.0 7r in steps of 0.1 tt. The maximum 
moves closer to zero energy with increasing 



to comment on a recent result obtained by Beri et al^ 
Using a different approach to calculate the conductance 
of a SC/HM point-contact, they find that generically 
G(eV = 0) = 0, at zero temperature. This agrees per- 
fectly with our results. One can show that below the 
gap, 



G(eV) oc \T*{e)\ 2 d v x{s,eV). 



(54) 



For T 



0, dyx = 2c5{e + eV) holds and since (1 — 
for e = we also find G(cV = 0) = 0. Note 
that ([53)) holds for arbitrary scattering matrices. Thus, 
this property is universal with the exception of $ = tt, 
ay = 0.57T where the denominator is zero for e = 0. 
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FIG. 13: (color online) The conductance G of an SC/FM 
point contact as function of contact voltage V (first row), 
and the corresponding SAR (second row) and AR (third row) 
contributions to it, for d = 0.1, 0.5, 1.0 Xf/^tt (left column) 
and d = 2.0, 3.0, 5.0 Af/27t (right column). The remaining 
parameters are E2 = 0.1i?F, E3 = 0.9 Ef, U+ = 1.1 Ef, 
U- = 1.9 E F , and a = 0.5 tt. 



B. Andreev conductance spectra 

As we have pointed out above, two competing Andreev 
processes participate in the presence of spin-flip scatter- 
ing, shown in Fig. [TT] Normal AR is suppressed as the 
polarization of the FM increases, since it requires one 
quasiparticle from each spin-band. SAR on the other 
hand takes two quasiparticles from the same band and 
thus dominates the spectrum for high polarization. We 
can define the corresponding contributions to the differ- 
ential conductance for each spin-band by 

^^dy^jte (F^3f^ 2 > + 
^f^=dv e -^lde(T^ 2 ri) + (55) 

and correspondingly (2 <-> 3) for band 3. The factor 
1/2 appears because the expressions in the integrand de- 
scribe only one of the two charges which are transferred 
by the process. The other charge is contained in X„ — x„ 
in Eq. (|44[) and cannot be disentangled from the one- 
quasiparticle transmission processes. The total contri- 



bution to the conductance is given by the sum over both 
bands: Gar = Gar,2 + Gar,3, Gsar = Gsar,2 + Gsar,3- 

In Fig. [T3] we discuss the results for the box potential 
using exactly the same parameters as in Fig. |TJ This cor- 
responds to a spin-polarization of the FM of P = 0.5. In 
Fig. [T3"k.b we plot the total differential conductance. For 
thin interfaces we obtain spectra with a rather conven- 
tional shape. The solid line in Fig. [T5k corresponds to 
a highly transparent interface, but still the conductance 
does not rise to a value close to twice the normal state 
conductance as in the conventional BTK picture. This is 
a direct result of the FM spin-polarization. Looking at 
the same line in Fig fT3T c) and Fig |T3T d). we see that the 
shape of Gar actually follows the usual trend, albeit with 
reduced magnitude, while Gsar gives almost no contri- 
bution in this case. The reduction of the Andreev con- 
ductance compared to the normal state conductance is in 
this case due to the spin-polarization of the FM. Only a 
fraction of the quasiparticles impinging the interface can 
undergo AR due to the reduced density of states in the 
minority band. 

As the thickness of the interface increases the conduc- 
tance contribution of SAR is enhanced and even dom- 
inates the sub-gap conductance for tunneling interfaces 
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0.8 0.9 1.0 1.1 1.2 0.8 0.9 1.0 1.1 1.2 



eV / A eV / A 

FIG. 14: The conductance G of an SC/FM point contact as 
function of contact voltage V , for (a) T = and (b) T — 
0.1 T c . In both cases, the values of E 3 = 0.2, ..,0.9, U- = 
E- +Ef are increasing in steps of 0.1 Ef from top to bottom. 
The remaining parameters are E2 = O.IEf, U+ = 1.1 Ef, 
a = 0.5 7T, d = 5.0 Xf/2-k. 



(Fig. [T3"H,D. This is because the magnitude of SAR is 
insensitive to the spin-polarization, as it takes two quasi- 
particles with the same spin from the FM. On the other 
hand it is very sensitive to spin-active scattering, which 
is why it is reduced for thin interfaces. We also see that 
as the transparency of the interface decreases, a sub- 
gap peak develops, as discussed in the previous section 
fFig. [T3"b) . However, the Andreev bound state stays close 
to the gap-edge in this scenario and smears out even for 
very small temperatures. 

In Fig. [Tla.b we plot the spectrum around the gap en- 
ergy for different polarizations, i.e. exchange fields, of 
the FM and a tunneling interface d = 5.0 Af/27t. Ap- 
parently, the sub-gap-peak moves to lower energies as 
the exchange field increases but also decreases in magni- 
tude. In any case, the peak is too small and too close 
to the gap-edge to be observable at finite temperatures 
fFigfT4b). This situation cannot be circumvented in the 
frame of the box-potential model, the reason being that 
one cannot obtain high mixing angles for reasonable pa- 
rameter ranges. Moreover, this situation is aggravated 
by the Fermi-surface average. As the mixing angle varies 
with the trajectory impact angle, the peak is broadened 
even at T = 0. This points again to the crucial impor- 
tance of the Fermi-surface geometry. If the Fermi-vector 
in the SC is considerably smaller than those of the FM 
bands, the scattering states which contribute to the cur- 
rent will be confined to a small range around perpendic- 
ular impact and hence a sharper peak structure can be 
expected. 

Finally, we show that even if this exotic feature in the 
conductance spectrum is not observable at finite tem- 
peratures, the impact of spin-active scattering can still 
be important. This holds in particular for FMs with 
high polarization, where SAR will naturally dominate the 
spectrum, if it is present. This can be seen in Fig. 1151 
where we plot the conductance for a highly polarized 
(P = 0.8) FM for a = 0.5 n and a = respectively. In 
the latter case, SAR cannot occur. If a = 0.5 it, the spec- 
trum is largely enhanced around the gap energy. This is 
not surprising, since SAR is mainly contributing in this 




FIG. 15: The conductance G of an SC/FM point contact as 
function of contact voltage V, for (a) T = 0, and (b) T = 
0.1 T c , shown for two values of a. The remaining parameters 
are E 2 = 0.1E F , E 3 = 0.99 E F , U+ = 1.1 E F ,U- = 1.99 E F , 
d = 1.0 Af/2tt. 



energy range. Even at finite temperatures an appreciable 
difference between the curves remains. 

Turning to the smooth scattering potential, we see that 
the situation changes fundamentally. We calculate the 
spectrum for the same set of parameters as in Fig. |TD] 
These results are shown in Fig. [T5] As we find a con- 
siderably enhanced mixing angle in this case, it is not 
surprising that the sub-gap peak is located far from the 
gap-edge if the potential is sufficiently smooth and may 
even be observed at finite temperatures. The width of 
this peak is directly related to the Fermi-surface average. 




eV/A 



FIG. 16: The differential conductance for the smooth po- 
tential model. The parameters are the same as in Fig. 1101 
The interface smoothness parameter a increases from back to 
front by steps of 0.1 Af. Temperatures are T — (top) and 
T = 0.1T C (bottom). 
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The calculations in Fig. [16] are for a tunneling limit situ- 
ation (t 2 < 0.01), and formula (f5"2"j) holds approximately. 
As one can see from Fig.fTOb. d sweeps through the whole 
range from to its maximum value as a function of the 
trajectory impact angle. This results in broadening and 
also implies that the Fermi-surface geometry may have an 
important impact on the shape of this bound state peak. 
For the particular geometry we consider here, with the 
Fermi surfaces of the FM-bands being both smaller than 
that of the SC, the mixing angle reaches for grazing 
impact. If however, the SC-band is smaller than at least 
one of the FM bands, this is no longer true as it can be 
seen in Fig. [4] Doing WKB calculations for different ge- 
ometries, we found that this may result in a kink at the 
tail of the peak, if ^ m ; n is large enough. 



C. Connection to the extended BTK-model 

The extension of the BTK-model to ferromagnetic 
point contacts proposed in Ref. [36| and further elabo- 
rated on in Ref. H3, was first used in Refs. I38ll39l to ex- 
tract the FM spin-polarization from the spectra of such 
contacts. Here, we show how this model can be ob- 
tained from our theory. The extended BTK-model char- 
acterizes interfacial scattering by a single parameter Z, 
which controls the transparency of the interface. Z is 
assumed to be independent of the transport channel. 
The spin-polarization of the FM is then taken into ac- 
count by noting that if P is finite, the transport channels 
can be divided into "non-magnetic" and "half-metallic" 
channel a 37 ' 39 , which is illustrated in Fig. [T] in this paper. 
This amounts to writing the conductance of the contact 
as a sum of the non-magnetic and half-metallic contribu- 
tion according to 3 ^: 



G = (1 - P C )G N + P C G 



H ■ 



(56) 



Here, the transport spin-polarization Pq was introduced: 



Pc 



N F2 vf2 - N F3 v F3 
N F2 V F 2 + N F3 v F3 ' 



(57) 



Gh is zero below the gap, since spin-flip scattering can- 
not occur in this model. This means that for |eV| < A 
one simply has the standard BTK formula reduced by a 
factor (1 — Pc)- The connection to our model is now es- 
tablished by making corresponding assumptions for the 
normal state scattering matrix of the interface. Since 
there is no spin-flip scattering, the matrix is necessar- 
ily diagonal, spin-mixing effects are obviously also disre- 
garded. Moreover, the fact that wavevector mismatches, 
let alone a spin-dependent interface potential, will intro- 
duce a spin-filtering effect is also not taken into account. 
Consequently, the whole scattering matrix is described 
by a single transmission parameter Tjy = t\ = t\. Eval- 
uating the corresponding expressions for T^' A is straight 



forward and yields: 

2-T^x n (e,eV)e(p F3 -\ Pll \) 



(58) 



1 



R 2 N -^(2e 2 -A 2 ) 



e < A 



with i?7v = 1 — Tn, and Q(p F3 — \p\\\) is the kinematic 
constraint for trajectories to be "non-magnetic". The 
total current density is: 



J(V) 



E 



eiV Fl? 



de(vr,j. 



(59) 



where we sum over the contributions of both bands for 
which an FS-average is calculated independently. 
Writing the FS-average explicitly, we get 



J = 



Note that 



de 



d Pp2 ~"F2j £ (PF2) _ 

(2^)3 \v F2 (p' F2 )\ + 

d 2 Pp 3 V F 3je(PF3) 

(2nhf \v FS {j/ F3 )\ 



Srj — 



«F,,,z>0 



d p Fn 



v Fri 



VF v (p' Frj )\ 



(60) 



(61) 



is exactly the area of the projection of the Fermi-surface 
onto the contact plane. Due to the kinematic constraint 
we have j(p F2 ) = j(PF3) = j(j>\\), if the parallel mo- 
mentum components of p F2 and p F3 are identical, and 
j(P\\) = f° r P|| > .PF3, which follows from equation (|58[) . 
This together implies that both integrals give the same 
contribution to the current which is not surprising, since 
Andrccv reflection induces the same current contribution 
in both bands. In the extended BTK-model case the cur- 
rent density j is not a function of pii , as T/v is not trajec- 
tory dependent. Assuming spherical Fermi-surfaces, we 
can hence calculate the FS average explicitly: 

-? eWF3^F3 f , . / , n _ rar ,s 

J =je z = / dej £ (e, V)e z . (62) 



The conductance is then given by: 



G% a 1 /Gn,o = Adyj /Gn.q — % 



N 



(63) 



where A is the contact area. Calculating T/v = 
J e de d v j E (e, V)/2e at T = yields the BTK formula^ 
with T N = 1/[1 + Z 2 } (note that we used dyx = 2e<5(er + 
cV) at T = at this point). G Nfi = (2e 2 Av F3 N F3 )/8 
is the contribution to the normal-state conductance of 
the non-magnetic trajectories. The corresponding term 
in (l56l) reads: 



(1 - P C )G N = 



2v F3 N] 



F3 



N F 2VF2 + N F3 V F3 



-G 



N- 



(64) 



To obtain the correct contribution to the normal-state 
conductance, Gjv must be related to the BTK-formula 
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by Gn = [e 2 A(Np2VF2 + -Wf3^F3)/8]7/v- Hence we have 
exactly 



N 



(1 — Pc)Gn- Analogously we can derive 
gquaM | e y| ^ ^ an( j recover the BTK result as well. 
We also obtain an expression for Gh for |eV| > A from 
our model by assuming a scattering matrix with n-f = 



T H e li} ' 2 : 
Gh/Gh,o = Th 



e""V 2 which implies T 3 = and T 2 = t\ 



(65) 



4T H (3 



2(/3 2 + !)-(£- lfT H -2cos?9Vl-T ff (/3 2 - 1) 



with j5 = eV/ v eV — A 2 . Comparison of this formula 
with that of Ref. H3 then shows that agreement requires 



costf = Z 



1 



l + Z 2 



1 



2(if/Z - 1) 



(66) 



(K — 2Z) 2 + 1, 

with T H = !/[! + Z 2 } and if is a parameter introduced 
in Ref. that we discuss in the following. For the 
sake of completeness, the corresponding contribution to 
the normal-state conductance is Gh.o = e 2 A(Nf2Vf2 — 
Nf3Vf3)/S. Apparently, a spin-mixing phase is manda- 
tory to reproduce the formula of Mazin et al. This result 
is not surprising, since the model used in Ref. [13 to calcu- 
late Gh necessarily introduces a spin- mixing effect, which 
is not true for the standard BTK model. The reason for 
this is that BTK assumes the same wavevectors in all 
channels and a non-spin-active interface. On the other 
hand Mazin et al. introduce different wavevectors by as- 
suming an evanescent mode in the minority band. This 
leads to the appearance of the quantity K = K,/k in their 
formula, where k controls the attenuation of the evanes- 
cent mode {e~ KZ ) and k is the component normal to the 
interface of the wavevector in the propagating channel. 
From our point of view this is nothing but a manifesta- 
tion of a spin-mixing phase, which is why we can only 
reach agreement by taking that into account. To make 
this point more convincing, we derived Eq. (|66[) from an 
explicit calculation of the normal state S-matrix using 
the same model as Refl37l. We match plane waves with 
wavevector k in all propagating channels and the same n 
as above for the evanescent mode of the minority band in 
the FM. The interface is modeled by a spin- independent 
delta-function with a weight factor. This yields the re- 
flection eigenvalues of the S-matrix on the SC-sidc n-j-, 
r U- By definition we have $ = arg[ri^r^] and find ex- 
actly Eq. (|66|) . In conclusion, we have shown here that 
earlier models for Andrcev reflection in clean ferromag- 
netic heterostructures are contained as limiting cases in 
our theory. As already noted in Ref. H3, the formula for 
Gh used by Soulen et al£2- was not obtained from a rig- 
orous calculation and is discontinuous at the gap-energy. 



polarized ferromagnets to study the conductance of 
SC/FM point contacts with a spin-active interface. We 
describe the interface by a microscopic model that ex- 
tends earlier models used in the description of Andreev 
reflection in such structures. Our main result are: (i) 
two types of Andreev reflection arise, one of them be- 
ing related to the creation of equal-spin triplet correla- 
tion. These processes depend differently on various prop- 
erties of the interface and bulk materials involved, (ii) 
the shape of the scattering potential has a pivotal im- 
pact in the magnitude of the spin-mixing effect. The 
usually assumed box-like or delta-function-like potential 
generically implies small mixing angles, (iii) we find spin- 
polarized Andrcev bound state peaks in the conductance 
of a point contact with a strong ferromagnet, that are 
more prominent for smooth interface potentials or a fi- 
nite magnetization near the interface in the superconduc- 
tor. The latter effect could be e.g. caused by the inverse 
proximity effect. Lastly, we would like to stress that the 
feature G(eV = 0) = for T = 0, which is universal for 
the spectra of half- metallic point contacts, may point to a 
criterion for identifying SAR in experiment at sufficiently 
low temperatures. 
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Appendix A: Scattering matrix parameters 

The general boundary conditions for the scattering 
problem in the quasiclassical theory of superconductivity 
have been derived in Ref. [ill These boundary conditions 
are formulated in terms of the normal state scattering 
matrix of the scattering region which has to be assumed, 
calculated from a microscopic model or fitted to exper- 
iment. In the case of a spin-active interface between a 
normal metal and a ferromagnet this matrix is a unitary 
4x4 matrix and one may ask for a set of parameters 
that describes the most general matrix uniquely and still 
allows for an interpretation of these parameters with re- 
spect to the physical problem in question. 



1. Singular value decomposition 



V. CONCLUSIONS 

In summary, we have used an extension of the qua- 
siclassical theory of superconductivity to strongly spin- 



Using a partial singular value decomposition (SVD) 
and the spectral theorem one can arrive at a decompo- 
sition of S that provides an appealing set of parameters. 
By partial we mean that a SVD is calculated for each 
block and not for the whole matrix, i.e. we have at the 
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outset: 



s = 



URV^ WT& 
WTZ^ -URV^ 



(Al) 



U, V, W, Z, W, Z, U and V are unitary and inde- 
pendent 2 x 2-matrices, while R, T, T, R are diagonal 
and contain the corresponding singular values. Such a 
decomposition is possible for any 4x4 matrix, which 
means that we did not exploit the unitarity of S so far. 
Exploiting unitarity we arrive at: 



S = 



U 
u 



R T 

T -R 



V* 

vt 



(A2) 



U, U, V, V are again unitary and independent. R and T 
contain the singular values of the composition and uni- 
tarity dictates i?i?t + TT^ = 1. To obtain a decomposi- 
tion which allows for a clcarcut interpretation in terms of 
scattering phases and spin-rotations, one has to continue 
decomposing U, U, V and V and eventually arrives at: 



S = 



Q 
Q 



$2 

l>i? 



Y 
Y 



y/l-TT 

-VI - TT 



Ft o 

ft 



$2 

$2 



Qt 
Qt 



(A3) 



This decomposition is written in terms of 2 x 2 blocks 
which are related to reflection and transmission, these 
blocks are matrices in spin-space. The central matrix 
contains the singular values of the partial singular value 
decomposition. These singular values relate to the trans- 
mission and reflection amplitudes of the interface but not 
in a simple way, since the outer matrices contain several 
rotations in spin-space. The outer matrices come in two 
flavors. The matrices Q, Q, Y, Y can be regarded as ro- 
tations of the quantization axis on either the left (Q, Y) 
or right (Q, Y) side of the interface. They have the 
structure: 



rot(a, ip) 



cos(a/2) sin(a/2)e iv 
sm{a/2)e- iip cos(a/2). 



(A4) 



The matrices $ 
complex phases 



$2 and ^2 are diagonal and contain 
(ijfs) 2 ). Their structure is: 



phased, ■&) 



e ir) e io- z ti/2 ^ 



(A5) 



Apparently r\ is a global phase and & a relative phase. 
The decomposition as it is presented here has 16 free pa- 
rameters which agrees with the maximum number of free 
parameters a unitary 4x4 matrix can have. However, 
we can now identify parameters which will be irrelevant 
for our problem. First we use the freedom of choosing 
an arbitrary quantization axis in the SC and put Q = 1. 
Secondly, we note that if the quantization axis of the in- 
terface does not rotate in the x-y-plane, we have S = S T 



and none of the rotation matrices defined above rotates 
in that plane. This implies Y' = Y T , Y' = Y T and 
Qt = Q T and also that <3> is real. We may hence absorb 
$ into T which amounts to having transmission eigen- 
values that can be negative. The transport properties, 
i.e. the current in this case, should also not depend on 
whether we extend the interface region arbitrarily far into 
the asymptotic region. This corresponds to the following 
transformation of the S-matrix: 



with 



S' = rjSrj 







P i(?72 +t73 + (m -V3 )fz ) /2 



(A6) 



(A7) 



Inspection of the boundary conditions shows in fact that 
both X v and T^' A are invariant under this transforma- 
tion. Considering this as another gauge transformation, 
one can eliminate the global phase in $ and use 772, 773 to 
obtain exactly the structure of (f3"2"]) in the transmission 
blocks. The reflection part on the SC-side reads: 



(A8) 



from which we conclude that the relative phase 
what is usually referred to as the spin-mixing-angle: 



is 



(A9) 



This is also the quantity which we plot in Fig. [4] [10] The 
necessity to have two additional mixing phases $2 and 
$3 comes about due to the additional rotations Q and 
Y. They are a function of all parameters which enter the 
transmission part. Thus a simple relation like (|A9|) does 
not exist in this case. 

The angle ay of which we make extensive use in the an- 
alytical discussion is associated to Y by Y = rot (ay, 0). 
In fact, the relations stated in (|3"4")l are given by YRY'. 

A fully analogous argumentation can be developed for 
the half-metallic case, however the corresponding scatter- 
ing matrix is 3 x 3 and thus all tilde-quantities are scalar, 
making them irrelevant. Furthermore, one can show that 
the ^> matrix is also just a scalar phase in this case and 
hence $2 fully accounts for the spin- mixing effect. So we 
have in the half-metallic case: 



# = 7^2=^$/2. 



(A10) 



This relation between and $2 is due to the fact that 
the evanescent solution in the ferromagnet is completely 
absorbed in the scattering matrix. 



2. Box potential 

For the special case of a box shaped potential we obtain 
analytical solutions for the scattering matrix, assuming 



17 



for the normal metal (superconductor in its normal state) 
a wave function of the form 



Sl- 



Ai+ \ -ik iz 
Ai- 1 



(All) 

with |fc||| 2 + k\ = 2mEp/h 2 , and in the barrier region 



where |fc||| 2 



k 2 



2m(E F - E 2 )/h 2 and |£||| 2 



2m(E F — E^)/h 2 , and in the case of two evanescent solu- 
tions, 



Jknf, 



D 2e -K2(*-a) 

D e -K 3 (z-a) 



(A15) 



(A12) 



with |feii | 2 — k± = —2m(U± — Ep)/h 2 , and with a certain 
spin rotation matrix U(a) that represents the misalign- 
ment of the barrier magnetic moment with the magne- 
tization direction in the ferromagnet by a misalignment 
angle a. The indices ± refer to spin-up and spin-down 
with respect to the misaligned spin quantization axis in 
the barrier. In the ferromagnet we can have, depending 
on the value of fc||, propagating or evanescent solutions 
in either of the two spin bands. In the case of two prop- 
agating solutions they arc 



V"2 

S3 . 



-iks (z—a) 



2_ p ik 2 (z — a) 



A3 iks (z — a) 



(A13) 



where |ft||| 2 + k\ = 2m(E F - E 2 )/h 2 and \k 



III 



2m(E F — Ej,)/h (in a more general model the masses 
on the two sides of the interface could also differ; we 
assumed them identical for definiteness). In the case of 
one propagating and one evanescent solution, 



S2 —ikf(z—a) 



A2 „ikf.(z—a) 

•/V2 



(A14) 



where |&||| 2 



2m(E F - E 2 )/h 2 and 



Nl 2 



2m(E F — Es)/H 2 . We then match the wave functions 
and their derivatives at z = (^tv arid ^b) and at z = a 
{^b and ^f), and eliminate the components D 2 and D3. 
The scattering matrix then is defined as the coefficient 
matrix in the relations 




= S 




(A16) 



for the case of two propagating solutions in the ferromag- 
net, 

A 1+ \ ( s 1+ \ 
Ai_ =5 si_ (A17) 

A 2 J V 32 J 

for the case of one propagating and one evanescent solu- 
tion, and 



Au 



S 



Sl- 
Sl- 



(A18) 



for the case of two evanescent solutions in the ferromag- 
net. 
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